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Abstract 

In this paper we propose some Hamiltonian characterizing the interaction of 
H ■ 

■ the two-level atom and both the single radiation mode and external field, which 



Q^i might be a generalization of that of Schon and Cirac ( quant-ph/0212068| ). We 



solve them in the strong coupling regime under some conditions (the rotating wave 
' approximation, resonance condition and etc), and obtain unitary transformations 



of four types to perform Quantum Computation. 

In this paper we consider a full model of the interaction of the two-level atom and 

both the single radiation mode and external field (periodic usually), which might be a 

generalization of that of Schon and Cirac 4J. We treat the external field as a classical 

one in this paper. As a general introduction to this topic in Quantum Optics see P], j2], 

0. Our model is deeply related to the (quantum computational) models 
*E-mail address : fujii@yokohama-cu.ac.jp 



(i) trapped ions with the Coulomb interaction, 

(ii) trapped ions with the photon interaction (Cavity QED). 

In our model we are especially interested in the strong coupling regime, [Oj, [IHI, 
The motivation is a recent interesting experiment, fTj. See jH] and as a general 
introduction^. 

In jO] and [12] we treated the strong coupling regime of the interaction model of the 
two-level atom and the single radiation mode, and have given some explicit solutions 
under the resonance conditions and rotating wave approximations. 

On the other hand we want to add some external field (like Laser one) to the above 
model which will make the model more realistic (for example in Quantum Computation). 
Therefore we propose the full model. 

We would like to solve our model in the strong coupling regime. Especially we want to 
show the existence of Rabi oscillations in this regime because the real purpose of a series 
of study ([IH], [H]; [IS]) is an apphcation to Quantum Computation (see [12] as a brief 
introduction to it). 

We can show the Rabi oscillations and obtain Rabi frequencies in this regime if the 
external field is constant. If it is not constant then the situation becomes extremely 
difficult. 

Let {(Ti, (72, cJs} be Pauli matrices and I2 a unit matrix : 



and CT+ = (l/2)(cri + i<Ji), a- = (1/2)((Ti — iai). Let W be the Walsh-Hadamard matrix 




(1) 



W 




P] is thick but strongly recommended 
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then we can diagonalize ai as cti = Wa^W ^ = ai = Wa^W by making use of this W . 
The eigenvalues of ai is {1,-1} with eigenvectors 

II I f I 1 I f I 1 

(3) 

Let us consider an atom with 2 energy levels Eq and Ei (of course Ei > Eq). Its 
Hamiltonian is in the diagonal form given as 



(4) 







This is rewritten as 



1/ ^ \ 
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where A = Ei — Eo is a energy difference. Since we usually take no interest in constant 
terms, we set 

Ho = -jas. (6) 

We consider an atom with two energy levels which interacts with external (periodic) 
field with gcos{u!Et)- In the following we set h = 1 for simplicity. The Hamiltonian in the 
dipole approximation is given by 

A 

H = Ho + gcos{ujt)(Ti = - yo-s + 9 cos{uEt)ai, (7) 

where o^^; is the frequency of the external field, g the coupling constant between the 
external field and the atom. We note that to solve this model without assuming the 
rotating wave approximation is not easy, see |H], JH], QH]; pH] - 

In the following we change the sign in the kinetic term, namely from —A/2 to A/2, to 
set the model for other models. However this is minor. 



Now we make a short review of the harmonic oscillator within our necessity. Let a{a^) 
be the annihilation (creation) operator of the harmonic oscillator. If we set N = a^a (: 
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number operator), then we have 

[A^, a^] = , [N, a] = -a , [a\ a] = -1 . (8) 

Let ?i be a Fock space generated by a and , and {\n) \ n G N U {0}} be its basis. The 
actions of a and on 7i are given by 

a\n) = ^/n\n — 1) , a^n) = \/n+ l\n + 1) , N\n) = n\n) (9) 

where |0) is a normahzed vacuum (a|0) = and (0|0) = 1). From state \n) for n > 1 
are given by 

|n) = ^|0) . (10) 



These states satisfy the orthogonahty and completeness conditions 

oo 

{m\n) = Smn , '^\n){n\ = l. (11) 

n=0 

Then the displacement (coherent) operator and coherent state are defined as 

D{z) = e^"^-""" ; \z) = D{z)\0) for z e C. (12) 

We consider the quantum theory of the interaction between an atom with two-energy 
levels and single radiation mode (a harmonic oscillator). The Hamiltonian in this case is 

i A + 
H = ujl2(S> a^a + —a^ ®l + gai® (a^ + a) (13) 

where u is the frequency of the radiation mode, g the coupling between the radiation field 
and the atom, see for example 0- 

Now it is very natural for us to include d?)) into (fT^ . so we present the following 

Unified Hamiltonian 

, , A 

H = ul2 ® a'a + giai ® {a^ + a) + —as 1 + g2COs{uj Et)cri ® 1. (14) 

Our Hamiltonian has two coupling constants. We note that our model is deeply related 
to the models 



(i) trapped ions with the Coulomb interaction 



(ii) trapped ions with the photon interaction (Cavity QED) 



|0) 



|1) 



This Hamiltonian is also related to the one presented recently by Schon and Cirac 0] 



P ( \ 

H = h co'ol2 ® a^a + gix) (a+ ® a + a- ^ a^) + 

v.m. V / 



2m 



-(T3 



(15) 



For the meaning of several constants see j4j. They have assumed the rotating wave 
approximation (see for example [Sj) and the resonance condition, and use a position- 
dependent coupling constant g{x), so their model is different from ours in these points. 
A comment is in order. Following |lj the Hamiltonian (fT^ might be modified to 

H = h 0^12 ® 0.^0, + 5'i(a;)cri ® (a"!" + a) + —0-3 (g) 1 + g2COs(iUEt)cri ® 1. (16) 

2m 2 

This model is a full generalization of (fTHjl . however we don't consider this situation in the 
paper. 

We have one question : Is the Hamiltonian ()14|1 realistic or meaningful ? The answer 
is of course yes. Let us show one example. We consider the (effective) Hamiltonian pre- 
sented by NIST group j3] , jH] which were used to construct the controlled NOT operation 
(see as an introduction). 



A 



(17) 



5 



We can show that under some unitary transformation the Hamiltonian ()17|) can be 
transformed to ()14|) with special couphng constants, [7j. This is important, so we review 
and modify [Zj. 

We set 2 A = iri{a^ + a) for simphcity, then 



(T+®e*^''("^+'') + (T_ 


^0 e2^ ^ 




^e-2^ ) 




^ e^ \ 




^ e-^ ; 




t/(77)((73® l)f/(r/)t 



i'q{a^ +a) 

A \ / 1 

1 



e 

e 

3-^ 



V 



1 1 



1 

Vo -1 



e^ 


1 1 
1 -1 



e^ 




(1^ 



where 



v/2 



e^ 




e^ — e^ 



V2 \ Q-A g-A 



(19) 



1 1 
1 -1 

and e"^ = D{irj/2) where -D(/3) is a displacement (coherent) operator defined by (fT^ . 
Then it is not difficult to show 

U{ri)^HU{ri) = ^^la ® 1 + CJ0I2 ® a^a+ ^c^i ® (-^a"^ + «a) +^(J3 (g) 1 - ^ffi ® 1. (20) 

To remove i in the term containing a we moreover operate the unitary one 



(I2 ® e*(-/2)^) UivYHUiv) (I2 ® e-*(^/2)^ 



= ^^l2 ® 1 + tuol2 ® a^a + -^ai ® (a^ + a) + c/aa ® 1 - — ai 1, (21) 

^ ^ 

where we have used the well-known formula 



see [T3]. Since U{rD can be written as U{vi) = (cr_|_ e^ + cr_ e '^) (W^ ® 1) , so if we 
write 

T{r]) = f/(77)(l2 ® e-^(-/2)A.) _ ® e^ + a_ ® e-^)(H^ ® g-'^'^/^)^), (22) 
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then we have 



T{ri)^HT{ri) = ® 1 + ujqU ® a^a + ^(^i ® {o) + a) + ga^®!- yai ® 1. 

Here we have no interest in the constant term, so we finally obtain 

H = T{7]) f^tuoU ® a^a + ® (a^ + a) + ® 1 - yai ® l| T{r])l (23) 

T{ri) is just the unitary transformation required. 

At this stage we would like to make a further generalization of the Hamiltonian ()14|) 
to make wide applications to Quantum Computation.® 

Let {K+, K_, K3} and { J+, J_, J3} be a set of generators of unitary representations 
of Lie algebras su{l, 1) and su{2). Then we can make the similar arguments done for the 
Heisenberg algebra {a\a,N}, namely (jH)) ~ (fT^ . see for example [TK] . 

We have considered the following three Hamiltonians in : 

(N) iJiv = Oa^a+ y(T3 O 1 + ® (a"^ + a), (24) 
(K) Hk = ujl2®K3 + ^a3®lK + 9(ri®{K+ + K_), (25) 
(J) Hj = iuU ® J3 + j(r-i ®lj + gai ® (J+ + J_). (26) 

To treat these three cases at the same time we set 



{L+,L_,L,} 



(N) {a\a,N}, 

(K) {K+,K.,Ks}, (27) 

(J) {J+,J-,J3} 



and 



Hl = UJI2 ®L3 + ®1l + gcri ® {L+ + L_). (28) 



Therefore the Hamiltonian that we are looking for is 
Full Hamiltonian 

El = UJI2 ®L-i + giai ® {L+ + L_) + ® 1l + g2Cos{ujEt)(Ti ® 11- (29) 



From now we would like to solve this Hamiltonian, especially in the strong coupling regime 

{gi > A). 

Let us transform into 



-f^L = I2 ® i^L^ + 0-1® {gi{L+ + L_) + g2COfi{ujEt)lL} + ^0-3 ® 1l 
A 



(30) 



The method to solve is almost identical to so we give only an outline. By making 
use of the Walsh-Hadamard matrix Q 

Ho = {W® 1l) [I2 ® ujLs + as ® {giiL+ + L_) + g2Cos{ujEt)lL}] (W'^ ® 1l) 
= E (|A)®e-^(^+-^-)){fiL3 + A^72CosKt)U}((A|®e^(^+-^-)) 

A=±l 

where |A) is the eigenvectors of ai defined in (jH)) and Q, x are given as 



(N) u, X = 2gi/uj, 

(K) Lj^l - {2gi/uy, a; = tanh-^(2(7i/c^), 

(J) LJ^l + (2(7i/^)2, x = tan-i(2(7i/^). 



(31) 



That is, we could diagonalize the Hamiltonian Hq. Its eigenvalues {En{t)} and eigenvec- 
tors {|{A,n})} are given respectively 



(E„(t), |{A,n})) 



(N) Q{-^ + n) + \g2Cos{uEt), |A) ® e'^^'^'-^V), 

(ir) r](ir + n) + A^2Cos(u;£;t), \X) ^ e-^^^+-^-^K,n) , 

(J) r](-J + n) + A5(2Cos(u;£;t), |A) (g) e^^^-^+^'^-^l J, ^) 



(32) 

for A = ±1 and n G N U {0}, where En(t) = En + Xg2Cos{LJEt)- Then Hq above can be 
written as 

Ho = T.T.Enm\n}){{X,n}\. 

A n 

Next we would like to solve the following Schrodinger equation : 



(33) 



To solve this equation we appeal to the method of constant variation. First let us solve 
= Hq"^, which general solution is given by \E'(t) = t/o(^)^0) where \I'o is a constant 
state and 

Uoit) = 5]5]e-^^*^"+^(^'^/"^)^'"("^*)^|{A,n})({A,n}|. (34) 

A n 

The method of constant variation goes as follows. Changing like \E'o — > ^o(^)) "we have 

ij^^o = jUo^as ® IlWo^o = ^Hf^o (35) 

after some algebra. We must solve this equation. Hp is 

= E E e^*(^--^")+^(^-'^)(^^/"-)^^"(--*) ({A, m} I (^3 ® 1l) \ {fi, n}) \{X, m}) ({/i, n} \ 
= EEe'^*"^"""^^'^^'"^"''^""^"''*^^((^N^"^^^"^~V)) l{A,m})({-A,n}|, (36) 

\ m,n 

where we have used (A|cr3 = (— A| and \n)) is respectively 



\n)) 



In the following we set for simplicity 



(iV) |n), 
{K) \K,n), 
(J) |J,n). 



e(t).,,^^. (37) 
Here we divide Hp into two parts Hp = Hp + Hp where 

H'f = j:Y.^"'"''^'H{n\e'^^'^-'-^\n)) |{A, n})({-A, n}|, (38) 

A n 

H"p = Y,Y. e'{*^(™-")+2^®W>((m|e^"(^+-^-V)) m})({-A, (39) 

\ m,n 

Noting ((n|e^*^^+~^"-'|n)) = {iji\e~^'^'"+^^'^\n)) by the results in section 3 of [13j, Hp can 
be written as 

H'p = Y. {He^^'-^-^-^n)) {e2^eW|{l,n})({-l,n}| +e-2^®W|{-l,n})({l,n}|} . 
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Here we want to solve the equation i{d/dt)'^Q = Hp^o completely , however it is not easy 
(see jH], [2n| and Appendix). Therefore we make a strong assumption. Namely we 
consider only the constant external field in (cu^; = in (jHTjl ). so 

m = g2t. (40) 

In this case we can solve the equation completely^. 

A comment is in order. For a short time span it may be not unrealistic to consider 
the above situation. Anyway, the author doesn't know whether it is reasonable or not. 

Therefore we have 

n 

Now it is easy to solve the equation i{d/dt)'^o = Hp^o, see Appendix. Next let us 
transform . 



Hp" = e'*^(™-") {((m|e"(^+-^-V))e''^'*|{l,W)({-l,^}| + 

((m|e-"(^+-^')|n))e'''^'*|{-l,^})({l,^}|}- (42) 



m,n 



For simplicity in the following we set 



then 





En,A - - 


((n|e"(^- 


^-^-)|n)), 






where 


K = X 


< (K) 


A n! /I 1 2y 

2 {2K).a ! 


-^-"F„( 


: 2K) where k = sinh(x 


. (J) 


2 p (1 




: 2J) where k, = sin(x) 



(43) 



(44) 



from the results in sectin 3.1 of |13j . 



^In the following the method to solve the equations is different from that of in which the 
Schrodinger cat states were used. However in this case we cannot use them, so the method becomes 
rather complicated. 
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Now let us solve (jH^ 

For that using the method of constant variation again we can set \l/o(^) as 

^o(i) = {iun,nan,i + M„,i2an-i)|{l,n}) + {Un,2ian,l + Un,22an-i)\{-l,n})} , (45) 



where from the appendix 
/ 

Unit) = 



\ Un,21 Un,22 









V 


exp < 






^ Q-2ig2t J 






(46) 



then we have a set of complicated equations with respect to {an^xit)}. However it is 
almost impossible to solve them, so let us make a daring assumption like JH] '■ for rn < n 

*o(^) ={{Um,liam,l + Mm,12am-l)|{l,"^}) + {Um,2ia'm,l + Mm,22am -l) | {" 1, "i}) } + 

{{Un,uan,l +Un,12an-l)\{l,n}) + (M„,2ia„,i + M„,22an,-l) | { " 1 , ?^}) } • (47) 

This ansatz is enough for our purpose because we are only interested in the Rabi oscilla- 
tions. Then after a long calculation we have 

T 

J. iy 


















Q-2ig2t 


1 







(4J 



where 



T^n = ((m|e^(^+-^-)|n)), T^n = ((m|e-^(^+-^- V)) 
and the equation with m < — > n in (jlSjl . We note Tmn = T^m and T^n = Tnm because x 
is real from (jHT|) . For the details of T^n or Tmn see P^ . 

We must calculate the right hand side of (jlH)). After some algebra 




exp < —it 



i?„^A 



T 



T 

J. rr 



-itXn 



^n' (49) 
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from pUj) and Appendix. Therefore the RHS of ()48|) is 



itQ{m—n) p 

2 



V 



T 

J. 



-itXn 



p-i 



where from Appendix 

Afc,+ = -92 + \/El^ + gl \k,- 

From now we divide into the four cases. 
Case (I) : (jKIH) can be written as 

^ it{n{m-n)+Xm,+ -X„ +} ^ 



-92 



El. + 9l 



I 



1 I . / Tnin 

m I _ 

ji(Am,--Am, + ) / \ T-i 

Here we set the resonance condition, namely 



\ 



-«t(An,- — A„,+ ) 



^{m - n) + Am,+ - An,+ = 0, 



or more exphcitly 



This is a rather comphcated equation. Then the above matrix becomes 



A, 



v 



r-1 

m 

„lt(Am,--Am,+ ) / \ 7-1 



-it(A„__-A„_+) 



By the way, since 



-'^X El A + 92 for k = m, n, 



(50) 



(51) 



we can take if (72 is large enough (so-called the rotating wave approximation^ which 
neglects fast oscillating terms). Therefore we have the time-independent (!) matrix 




T 



I 



\ 



r 



-1 



■^in many papers in Quantum Optics this assumption has been used without showing that the frequency 
in the model is large enough. However it is not correct. 
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2 " 



7« 



r 



-1 



where 



In 



A 



-T 



n,A 



m, A 



+ -^m,A Y '^n,+ + -^n,A \/ -^"1,+ + -^m.A \/ '^n,+ + -^n,A 

1 1 



-T 



A 



y2 



-^m,A 



from = T„m. Therefore 7mn = Inm = 7 G R. As a result we obtain the equations 



(52) 



'~dt 



.d_ 



\ O'm-l 

0-11,1 
\ On-l 







Om,l 
\ Om-l 



(53) 



or if we introduce the compact notation = (am.,i, 0^,-1)* then we have the clear 
matrix equation 



. d ( 



A 
~2 



I 



V 



/ 



O2 



(54) 



Therefore the solution that we are looking for is 
a„(t) 



, . A 

exp < —it— 



I 



/ 



V7^ 



am(0) 



a„(0) 



• (55) 



O2 7^ 2 

02 y J V 

The Rabi frequency is just A7 (compare with JHDj which is also rather complicated. 

Next we consider the remaining three cases : The arguments are almost identical, so 
we only give the results (we leave them to the readers). 



Case (II) : The resonance condition is 



VL{m - n) + Am ,- - A„ _ = 0. 



(56) 
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The solution that we are looking for is 
a„(t) 



A / O2 7^ 



exp { -it— 



am(0) 
a„(0) 



where 7 is 



7 



The Rabi frequency is A7 (compare with |.13J. 



Case (III) : The resonance condition is 



n(m - n) + Am,+ - A„ _ = 0. 



The solution that we are looking for is 
a„(t) 



/ 



exp < —it- 



V 



A ( O2 7cr+ 

7Cr_ O2 



-1 



am(0) 
a„(0) 



where 7 is 



7 



(-E'n,A-^nmAm,+ ~l~ -S'm,A-^m7iAr 



The Rabi frequency is A7 (compare with [T^). 



Case (IV) : The resonance condition is 



f2(m — n) + Am,- — A„,+ = 0. 



The solution that we are looking for is 



/ 



V a„(t) 
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exp < —it 




-1 



am{0) 
V a„(0) 



(63) 



where 7 is 



7 



(64) 



The Rabi frequency is A7 (compare with ^13]). 



On the ansatz (jlTj) we solved the Schrodinger equation in the strong couphng 
regime (!) under the resonance conditions and rotating wave approximations, and ob- 
tained the unitary transformations of four types which are a generahzation of jJSl ■ They 
will play a crucial role in Quantum Computation. 

On the other hand we in this paper solved the special case ()4()j) . however we would 
like to study the general case (|H7j). At the present it is almost impossible (see Appendix). 
We will treat this case in a forthcoming paper. 

By the way, we considered one atom with two-level, so we would like to generalize 
our method to n atoms (with two-level) interacting both the single radiation mode and 
external periodic fields like {n atoms trapped in a cavity) 



|0) |0) 



|1) |1) 



|0) 



Then the Hamiltonian may be 

n ^ n n 

HnL = ojIm 0Ls + giJ2 ^1^'^ ® (^+ + ^-) + ^ ^ ^^^^ ® 1l + ^2 E cos{ujjt + 0j)(^P'^ ® 1l, 

(65) 

where M = 2" and ai^^ {k = 1, 3) is 



I2 ® (g) I2 



I2 (j — position). 
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In the near future we will attempt an attack to this model. However according to increase 
of the number of atoms (we are expecting at least n = 100 in the realistic quantum 
computation) we meet a very severe problem called Decoherence, see for example fOj and 
its references. The author doesn't know how to control this. 

A generalization of the model to N-level system (see for example ^3], ^Tj, ^H], [23]) 
is now under consideration and will be published in a separate paper"^. 

Acknowledgment. The author wishes to thank Yoshinori Machida for his warm hospi- 
tality during the 11-th Numazu Meeting held at Numazu College of Technology (6- 
8/March/2003). 



Appendix : Some Useful Formulas 



In this appendix we solve the following equation 



i— V' = aHip, 
at 



where a is a constant and 



H = e^'^'\l){-l\+e-^'^'\-l){l\ and ^ = a{t)\l) + b{t)\-l) . 



Then it is easy to get a matrix equation on {a, b} 



d 
I — 
dt 



/ 






a 










V 







e 

-2i0t g 



2iet 



The solution is easily obtained to become 



V ^0 



(66) 



(67) 



(69) 



"'The author believes that it is important for us to consider the N-level system to prevent the deco- 
herence problem 



16 



where (ao, 60)* is a constant vector and 

m 

If we set 

m 



fl 


\ 


exp. 


1- 


^ a \ 


1 


\ 


„-2iet 

^ J 






-26 ) 





d 

i—U = HU. 
dt 



(70) 



Uu U12 



\ 



(71) 



U21 U22 I 

(2 corresponds to —1) then above can be written as 

i) = (i^iiao + uviho) 1 1) + (ii2iao + ^^22^0) I -1) (72) 
with constants {ao, &o}- 

In the method of constant variation in the text we change hke ao — > ao(0 and 60 — 
6o(t). 

Let us make some comments. For 

\ 



A 



we can easily diagonahze A as follows : 



a 
a -26 



(73) 



A+ A- 



A. 



A_ 



A+ A- 



(74) 



where 



A 



+ 



-6 + Vd^Tc^, X- 



Therefore we obtain 



V{t) = e-''^ = 



A+ A_ 



-itX- 



A+ 



(75) 



We note the formula in the text. 



A+ A- 



A+ A- 



A_ —a 



A_ - A 



+ \ -a -A 



+ / 



(76) 
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For the simplicity we set 



( ^ 



In the text this is used as 



( Ek,A 



e 0(2). 



(77) 



(7J 



with 



A, 



-92 + \Jgl + El^, h,- = -g2-\jgl + 

for k = m, n. 

Last we make a comment. The fact is that we wanted to solve the following matrix- 
equation instead of (jUHj) 

/ . 







/ 





g2ie(t) 




= a 












V 


g-2ie(t) 






where 



We here note 



. d 

'dt 



(79) 



UJ 



6t as uj 0). 



g2ie{t) ^ ^ J„(0/cj)e"'' 



^0) 



neN 



where J„(x) are the Bessel functions. 

However we cannot solve this equation completely, see jH], [H], [2D1, [HI- This is just 
the bottleneck. 
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